d-Wave Pairing State in Terms of the Zhang-Rice Singlets 

Akinori Tanaka0 

Department of General Education, Ariake National College of Technology, Omuta, 

Fukuoka 836-8585, Japan 

(April 25, 2007) 



Abstract 

In cuprate superconductors doping is believed to create holes on the O-sites, which 
couple antiferromagnetically with holes on the Cu-sites to form the so-called Zhang- 
Rice singlets. Here we study a d-wave pairing state based on the Zhang-Rice singlet 
states. Upper and lower bounds of an off-diagonal long-range order parameter with 
d-wave symmetry for this state are estimated. We also introduce a concrete model with 
on-site Coulomb repulsion and kinds of antiferromagnetic interactions whose ground 
state is this <i-wave pairing state. 
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1 Introduction 



The mechanism of high-T c cuprate superconductivity has been attracting much interest since 
it is discovered in 1986 p. In cuprate superconductors, electrons (or holes) in the CuC>2 
planes play major roles, and the importance of the Coulomb repulsion at the Cu-sites is 
emphasized from the beginning [2j[3jll]. However, theoretical understanding of its effects on 
the superconductivity is still limited and is being a challenging problem in condensed matter 
physics. 

Most theories which start with viewing cuprate superconductors as doped Mott insulators 
are based on the so-called Zhang- Rice singlet states j5] . In the undoped case, where there is 
one hole per Cu-site in Cu02 planes, the cuprates exhibit insulating antiferromagnetism due 
to the strong Coulomb repulsion at the Cu-sites. When the system is doped, additional holes 
are created on the O-sites. Because of a superexchange antiferromagnetic interaction, each 
of the holes occupies a quasi-localized state on the four nearest neighbour O-sites around a 
Cu-site, forming a local spin-singlet with the hole on the central Cu-site. This singlet is now 
referred to as Zhang- Rice singlet. The Zhang- Rice singlets become charge carriers moving 
through the Cu0 2 plane and condense into a superconducting state. 

This scenario is usually examined by using the t-J model which is a single-band effective 
Hamiltonian with antiferromagnetic interactions between nearest neighbour holes on the 
Cu-sites [5|. Despite its simple form, however, it is a formidably difficult task to rigorously 
analyze the t-J model, and whether the model really describes the cuprate superconductivity 
has not yet been clarified. In the current situation, we think that a rigorous establishment of 
occurrence of a superconducting state based on the Zhang-Rice singlets in a model with the 
Coulomb repulsion and antiferromagnetic interactions, even if it is apart from the t-J model, 
certainly gives us an important step toward understanding of the cuprates superconductivity. 

In this paper, we study a simple <i-wave pairing state expanded in terms of the Zhang- 
Rice singlet states. It is shown that the pairing state is regarded as a condensed state of the 
Zhang-Rice singlets in the background of a resonating-valence-bond state consisting of holes 
at the Cu-sites. We estimate an upper bound on an off-diagonal long-range order (ODLRO) 
parameter with ci-wave symmetry for the pairing state as a function of doping concentration 
< 8 < 1. It is found that an upper bound has a dome structure with a maximum at 
5 = 0.5 and becomes zero at 5 = 0, 1. We also estimate a lower bound on the ODLRO 
parameter and show that ODLRO exists for sufficiently large doping concentrations. We 
then introduce a model with on-site repulsion and kinds of antiferromagnetic interactions, 
and show that the pairing state is a ground state of this model. A related model with 
infinitely large on-site repulsion at the Cu-site is analyzed in Ref. [6]. This model, however, 
has following disadvantages: the Hamiltonian does not have spin rotational symmetry, and 
its exact pairing ground state has less relevance to the Zhang-Rice singlets. Although the 
present model has still somewhat artificial aspects, it is for the first time that the pairing 
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state with <i-wave symmetry which is written explicitly in terms of the Zhang- Rice singlet 
states is realized as a ground state of the concrete Hamiltonian. 

This paper is organized as follows. In the next section we prepare some notation and give 
a definition of the Zhang- Rice singlet states. In section [3l we introduce a two-electron state 
with d-wave symmetry, and, on the basis of the Zhang-Rice singlet states, we construct a 
pairing state in which many electrons condense into this two-electron state. In section HI we 
discuss an expectation value of an order parameter with <i-wave symmetry for the pairing 
state. An upper bound for the order parameter is obtained in this section and a lower 
bound, whose estimation needs somewhat technical calculations, is obtained in section [6j 
In section |5] we introduce a Hamiltonian whose ground state is the pairing state which we 
construct. In the final section, a summary and some remarks are given. In Appendix lAl we 
show that the pairing state is non-vanishing. 

2 Zhang-Rice singlet states 

We start with the definition of a lattice. With even integers L\ and L 2 , let 

D = ([l,L 1 ]x[l,L 2 ])nZ 2 , (2.1) 

which represents a collection of the Cu-sites. Let 5 1 = (1, 0) and 5 2 = (0, 1). We define 

p = { u | u = x + 5 l /2, I = 1,2, x G D}, (2.2) 

which is the collection of the mid-points of the nearest neighbour bonds in D and corresponds 
to the O-sites. Then we consider the lattice A = D U P, which mimics the CuC>2 plane. (See 
Fig. [TJ) For a technical reason we impose periodic boundary conditions on A. For later use, 
we introduce further the following sublattices of D: 

D Q = {x | x = (x\, x 2 ) G D with x\ + x 2 being odd}, (2.3) 
D c = {x | x = (xi, x 2 ) G D with X\ + x 2 being even}. (2.4) 

Next we introduce fermion operators which annihilate or create holes with spin o =f , \ 
at sites in A. Any states with the number A r h of holes can be constructed by operating these 
operators on a state $o with no holes on A. By d x ^[d\ a ) and p u ,a{Pu o-)> we denote the 
annihilation(creation) operators of holes at x G D and u G P, respectively. As mentioned in 
section 1, each hole additionally induced in a CuC>2 plane with 1 hole per Cu is considered 
to localize well at the four nearest O-sites of a Cu-site because of the antiferromagnetic 
superexchange interactions between Cu- and O-sites. To describe this localized state on the 
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Figure 1: The lattice structure. The solid and open circles indicate the Cu- and O-sites, 
respectively. 

O-sites we introduce the following operators for each x G D [7j: 

fx,a — ~Z ^ Pu,a- (2-5) 

ueP;\u-x\=l/2 

As is easily seen, the annihilation operator f Xj(T and the creation operator f\, defined by 
(I2.5p do not anticommute when \x — x'\ = 1, implying that the single-electron states cor- 
responding to (12. 5j) are not orthogonal. To avoid technical complexities arising from this 
fact, we consider corresponding Wannier states. To do so, we introduce the fermion operator 
/* = (1/vZ?) J2x£D elnS x Px+s 1 /2,a and the reciprocal lattice 

K = {(-jr-ni, ^-n 2 ) | meZ, - L t /2 <n t < L t /2 with I = 1, 2} , (2.6) 

and then define / fej(7 = (l/y/\D\) J2x€D fx^~' lkx for k G /C\{(7r,7r)} and = We 

normalize the /-operators as afc jCT = fk,a/ \\ fk II? where the normalization factors are given 
by 

{1 if k = (tt, it), 

^1 + 1 (cos ki + cos k 2 ) otherwise. 

The fermion operators corresponding to the Wannier states are defined by 

= -i= J] «fc, CT e ife - x . (2.8) 
fce£ 

The a-operators defined as above approximate the /-operators well, and satisfy the canonical 
fermion anticommutation relations {a xa) a yT \ = { a x,a, a>y,r} = and {a' xa , a VtT } = 5 a}T 5 x ^ y 
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Figure 2: The solid and open arrows indicate spins of the holes on the Cu- and O-sites, 
respectively, (a) In the case of iVh = \D\, every Cu-site is occupied by one hole, (b) When 
iVh is greater than \D\, every Cu-site remains to be occupied by one hole, and additional 
holes are created on the O-sites. Each hole on the O-sites occupies a quasi-localized state, 
which is indicated by dot lines, and couples to the hole at the central Cu-site to form the 
Zhang- Rice singlet. 



for cr, t =|, | and x,y G D. In the rest of this paper, we consider the Zhang-Rice singlets by 
using the a-operators, instead of the /-operators. 

The Zhang-Rice singlet around a Cu-site x is formed by holes occupying a\ a and d) XT . 
This singlet is represented by the two-hole creation operator 

^ = < T < 4 + < T 4 4 - (2-9) 

We assume that, in the case where the hole number is \D\, each hole occupies a Cu-site. 
Any | D | -hole state is then expressed by a linear combination of YixeD <4 a x ^o with a x =f , I 
(Fig. 12(a)). We furthermore assume that N holes added in this state form Zhang- Rice 
singlets. Then a (\D\ + iV)-hole state with < iV < \D\ is written by using a set of states 

fn<J f n (2.10) 

\x&A / \y<=D\A ) 

where A is a subset of D with \A\ = \D\ — N and its compliment D\A is a collection of sites 
where the Zhang- Rice singlets are formed (Fig. [2(b)). Noting the relation 

4, CT $o = -sgn[a]a Xi _ ff ^$o (2.11) 
where sgn[a] = + if a =f and sgn[cr] = — if o =|, we find that f)2.10p is rewritten as 

(n«^Po = ^4,<x A , (2.12) 

\x£A ) 
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with 

*o=[II^Wo (2.13) 

up to a sign factor. Here a a is a short hand for a spin configuration (cr x ) xe A- We write Sa 
for the collection of spin configurations {(<j x ) xG a I o~x =T> I? x e ft is easy to see that 
(^A,cr A , ^b,t b ) — 2 |D| ~ |j4| xL4 = B]x[cr A = t b ], where x[event] = 1 if 'event' is true and 
otherwise. Thus the collection of states 

{*a,<t a I A CD, <r A eS A } (2.14) 

is orthogonal. For \D\ < iV h < 2\D\, let H^ s be the Hilbert space spanned by the basis 
states {^a,ct a } with \A\ = 2\D\ — N^. The Zhang-Rice singlet states are defined to be states 

in h zrs- 



3 rf-Wave Pairing State 

Assuming that takes an even number in \D\ < < 2\D\, we consider a <i-wave pairing 
state in the Hilbert space H^ s . Let us define a pair operator £ by 

C= E (cosfci - cosA;2)a_ fcii a fciT . (3.1) 

k=(ki,k 2 )eK 

Recall that dk, a = i)-/ \J\D\) J2 x ,eD a>x,cr&~ lk ' x are the Fourier transforms of a XjCT . This operator 
creates an electron pair with d-wave symmetry. For the hole number iVh = \D\ + N with a 
positive even integer N, one of the simplest d-wave pairing states in H^ s is given by 

$ p = (Cf P *o (3.2) 

with the number of pairs N p = (\D\ — N)/2. Here we note that 

a Xti a x ^ = (3.3) 

since [a X) ia x ^, ipl] = —al +a x j — d* x ,a X) i, so that $ p is actually expanded in terms of (I2.12p . 
In Appendix |A] we show that $ p is non- vanishing. 

In order to see the real space representations of ( and <£> p , we define the following opera- 
tors: for x G D 

b x ,a 

and for x,y E D 



y£D;\x-y\=l 



(3.4) 



~-e m5 ix y \a xd a yA + a y>l a xA ) 



(3.5) 
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The operator 0?^ y y corresponds to a two-electron singlet state formed by electrons on the 
O-sites around Cu-sites x and y. By using these operators we can write ( as 



or 



C — E °aa^,T — E ^B,l a *>T 

xeD x£D 

c= E 

{a;,j/}eB 



(3.6) 



(3.7) 



where 



B = {{x,y}\x,yeD, \x - y\ = 1} (3.8) 

is the collection of bonds in .D (we assume that {x, y} = {y, x}). Let C(B) be the collection 
of subsets B of B such that no two elements in B share the same site. Substituting (13.71) 
into (13.21) . and noting the relation (13. 3p . we obtain 



$™ = N, 



e n - 

BeC(B);\B\=N p {x,y}£B 



(3.9) 



Therefore, the pairing state $ p is regarded as a nearest-neighbour resonating-valence-bond 
state (which is a linear combination of products of two-electron singlets) consisting of elec- 
trons on O-sites with the background of the fully-filled Zhang-Rice singlets. 

Finally let us see the form of $ p in terms of the hole creation operators. Let n d x a = d\ a d x ^ a 
and define Vd — Y[ x &d{^- ~ n x,-\ n x,\)i which is the projection operator onto the space without 
double occupancies of holes at the Cu-sites. By using this projection operator we can rewrite 

as 

*o = tL'Pd f E Vl) ' ' $o = t^-Pd ( E + alAk,l)) $ o, (3.10) 

XxeD ) \ u \- \keK j 

where d k a = (1/v/ \D\) J2xeD d X:(T e~ lk ' x . Then, noting two commutation relations 



"-/.-. I «'k,T> E (4>T a -P,J. + °i,l$-P,l 
pe/c 



T a fcjT + d_ k ^ 



a -k,i 



(3.11) 



and 



' I E V"'' 

peK. 



(3.12) 
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we obtain 



$ p = ~m V ° E ( cos k<2 ~ cos fc i)4,T^fc,J, 





E II (^}) t ) $o (3.13) 

BeC(e)B;|B|=Ar p {x,i,}sB 



with 



^U}) ] = ^ sl<x ' y) (dl,A,i + 4,t4J- ( 3 - 14 ) 

To get the second line in (13.131) we used 7- > d4,t4 ) | = 0- The operator (^j^yi)^ corresponds 
to a two-hole singlet state formed by holes at the Cu-sites. From expression (13.131) we 
find that the state $ p can be regarded also as a projected state in which the Zhang-Rice 
singlets condense and the remaining holes at the Cu-sites are forming nearest-neighbour 
singlet states. 

Here it should be noted that, despite of the form (I3.13p . $ p does not exhibit long-range 
order associated with the Zhang- Rice singlets. In fact, it is easy to see that 

$ p ,^$ p ) = (3.15) 

for x ^ y, since there is no charge fluctuation on the Cu-sites. It is important to consider 
long-range order associated with movable holes on the O-sites, which we discuss in the next 
section. 



4 Order Parameter 



In this section we estimate the value of a ci-wave order parameter for the state (13.21) . Let 

1 ^ ,„ 1 



A 



\D 



E r {xM 

{x,y}£B 



\D 



<■ 



(4.1) 



We then define 



($ p ,AtA$ p ) 

\ <*p,*p> : 



lim fi AN , 

|D|,JV->oo 
N/\D\=S 



(4.2) 
(4.3) 



where the limit is taken with N/\D\ kept fixed to 5. This order parameter measures a long 
range correlation between spin-singlet pairs corresponding 4>{ x ,y}- 



We firstly show that 



$ p , AtA$ p ) = ^f±± £ E ( $ p> 4,^,.^,-A,-^p) (4.4) 



which is crucial for our estimation of fis (recall N p = (\D\ — N)/2). To see this, we observe 
that 

((c)^o,4A(cf p+1 *o) = ((o N <i n ^o^i^o^ 1 ^) 

\ yeD\{x} I 

= (Ap+i)((oMi n ^0,6^(0^4^ n v>j* ) 

= (A^ p + l)((C) JV ^o,< T ^, T ^,|^i(C) iVp ^o). (4.5) 
To get the second line we used the commutation relation 

[C4,«r> a l,M = sgnM^ (4.6) 

which immediately follows from the real space representation (13.61) of (. Then, by using the 
spin-rotation symmetry for <3> p , (I4.4p follows from (14. ip and (14.51) . 
By noting the inequality 



< (4-7) 



and the fact that the number of a-holes (which are holes in the state corresponding to the 
a-operators) is exactly N for <3> p , we find that //a, at is bounded from above as 



fJ'A.N 



< 



/iVp + l\ (N_ 
\\ \D\ )\2\D\ 



The limit is thus bounded from above as 



fis < -^5(1 -8). (4.9) 

As for a lower bound for fis we have the following result. Let | < S < 1. Then we have 
that 

H6>y-ysI(6)(l-6), (4.10) 
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where 75 = 2(8(5 _ 7 ) and 

J W = 7TY2 [ , ^) X [e b (k) < e s ]dk (4.11) 

with 6b(k) = (cos/ci — COS/C2) 2 . Here es is determined by 

6 = 7tW, Xh(fc) < e 5 ]d/c. (4.12) 

The inequality (I4.10p means that the state $ p exhibits ODLRO with d-wave symmetry for 
I < 5 < 1. The calculation for this bound is somewhat complicated and technical. We defer 
the proof to section [61 It should be noted that the above lower bound is not optimal at all 
and never means that there is no <i-wave pairing order in a low density region of holes. It is 
desirable to obtain an improved bound in the future. 



5 Hamiltonian with Ground State & p 

So far we have constructed the pairing state <3> p with <i-wave symmetry and studied its 
properties. In this section we propose a Hamiltonian if on A whose ground state is given 
by $ p . 

Let us define the number operator n^ a , with a =f , I, and the spin operators S* a , with 
a = 1, 2, 3, corresponding to the a-operators by 

1 

S x,l = 2^.T fl ^i + c 4,j. a ^t); ( 5 - 2 ) 
■y 

Sx,2 — 7^( a x ,r a x,i - a x,i a x,i), (5-3) 



2i 



^(4,T a ^T - 4,i a x,i)- ( 5 - 4 ) 



We also define 



The number and the spin operators for the b- and the (i-operators are defined similarly. By 
using these operators, the Hamiltonian H is defined as follows: 

H = H + H 1 (5.6) 
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with 



Ho = -e d J2n d x + Uj2n d x y Xtl + J J2S a x .S d x , (5.7) 

x£D x£ x&D 

3 

+J X £ {SI -S d x + S x .S d x + S a x .S b x - -nl ■ n b x ) . (5.8) 



x<=D 



Here, all the parameters, £ d , U, Jo and J±, are positive, and e d is assumed to take values 
in | J < e d < | J + U. It should be noted that one can rewrite H by using the d- and the 
p-operators, although it has a somewhat complicated form. It is also noted that we do not 
take any peculiar limit, such as U — > oo and Jo — > oo, and thus /J acts on a whole Hilbert 
space constructed by the d- and the p-operators. 

We shall show that the lowest energy of H for the hole number N h = \D\ + N with 
< N < \D\ is e = —e d \D\ — §JoiV, which is attained by the Zhang- Rice singlet states in 

n ZRS- 

Let N d be the eigenvalue of J2 x &D n t, the number of d-holes. Since N£ is a conserved 
quantity for H , it is convenient to decompose the iV h -hole Hilbert space into the subspaces 
with fixed N d . We denoted by the subspace with fixed N d and by E(N d ) the lowest 

h 

energy of Hq for the states in H^J. 

h 

Let us examine each term in Hq. The eigenvalue of the first sum in Hq is —e d N d for the 
states in H^. The lowest eigenvalue for the second sum is zero which is attained by the 

states without doubly occupied d-states. The eigenvalues of JoS® ■ S d are — | J , and | J . 
We have eigenvalue — |J when each of the d-state and the a-state at site x is occupied by 
one hole and furthermore the two holes in these states form the spin-singlet state. 

It immediately follows from the above observation that E(\D\) = eq, which is attained by 
the states in H^ s C H^. In the case < Ng < \D\, noting that there are JVJ = \D\+N-N£ 
holes on the O-sites, we have 

E(K) = s d N d -p mm(N d ,K) 

> e + 3 -J {K-min(NtK)} 

> e . (5.9) 

Here the second line follows from the assumptions < |J < e d and N d < |-D|(or N < jV£), 
and the third line follows from N£ > min(A^, N%). In the case \D\ < < \D\ + N, noting 
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that there are, at least, (N£ — \D\) doubly occupied d-states, we have 

E(K) = -e d N* + U(N*-\D\)-h NV 
= e +(h + U-e d )(N*-\D\) 



> e . (5.10) 

Here the final inequality follows from the assumptions Ed < + U and \D\ < N^. As a 
result, we have E(N^) > e for ^ \D\, which proves the claim. 

We have shown that the lowest-energy states of H Q are the Zhang-Rice singlet states in 
H^ s . In the following we shall show that Hi is positive semi-definite and $ p in H^ s is its 
zero energy state. This implies H = H + Hi > e and H$ p = e $ p . We thus conclude that 
$ p is a ground state of H . 

By a straightforward but somewhat lengthy calculation, one finds that H x is rewritten 
as follows: 

^i = ^iEEE + K™i(K™tf] (5-11) 



° xeDm=ll=l 



with 



and 



K l,i = b U a ^i d ^ 

K l,2 = ( 6 l,T a *a d *,T + 4,T a ^.T^,4 - b x,{ a x,l d x,l) > 

K l,3 = ( 6 l,T a *.T d *,T ~ b l,l ax >l d x,1 ~ b l,i a xA d x,i) ' 

k I,a = ~ b l,i a x,ld x , h 

K l,i = 4,T 6 *.l d x,i» 

K l,2 = ( a l,i b x,i d x,i + a l,i b x,id xd - a\ A b xA d x ^ , 

K x,3 = ( a l,i b x,i d x,i - 4,i 6 *.i d *.T ~ 4,i 6 *.T d aa) » 

^,4 = -4,i 6 *.T d x,T- 



5.12) 
5.13) 

5.14) 
5.15) 

5.16) 
5.17) 

5.18) 
5.19) 



It follows from this representation that Hi is positive semi-definite. Therefore, the lowest 
energy of Hi is greater than or equal to zero, and any zero energy state $ of Hi, if it exists, 
must satisfy (Kffl® = and K™jfi = for all m = 1, 2, I = 1, . . . , 4 and x e D. We shall 
prove that $ p indeed satisfies these conditions. 
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We start with the case of m = 1 and 1 = 1. By using the commutation relation (14.61) . we 
have 

(/V,'.:K = 4,i<; 6 *,tC = <i( 6 *,T) 2 + C<i4At = U/V,',:) • (5.20) 
This together with (K* J VI = 0, which follows from (aJ. CT ) 2 = (d) xa ) 2 = 0, leads to 

(Cf P *o = (Cf P «i)^o = 0. (5.21) 
From a x ^d x ^ipl = — da.^4 14 + a *,i4 t4 we immediately obtain 

(Cf P *o = 6l, t (C)" p = 0. (5.22) 

We thus conclude (K^y^p = K\ x $ p = for all x in D. 

Let us consider the cases of m = 1 and I = 2,3,4. Define spin-lowering and raising 
operators as 



S~ = E(4^,t + <At)> (5-23) 

x&D 

S + = E(4,T«-l + 4, T <l)- (5-24) 



x&D 



From the results for I = 1 we have 5 + (JiT* jjt^p = 0. It is easy to see that S + (K^ oi y = 

Vs(Ki 2 y + (K i x>1 ys + , s+(Ki i2 y = 2(k^ + (Ki >2 ys+, and s+(Ki 3 y = Vs(Ki,y + 

(K] c 3 yS + . Substituting the first relation into S + (K] ci y^ p = and noting S + § p = 0, we 
find (K^y^p — 0. Repeating the same argument, we have (iCj 3 )^ p = (fCj 4 )T$ p = 0. By 
using S-K^ p = 0, S-Kj^ = -yftK^ + K^S", S~K^ = -2i^ 3 + K^S", S'K^ = 
-V3Kl A + K^ 3 S~, and S~<& p = 0, we similarly obtain K^ p = for I = 2, 3, 4. 

Proceeding in the same way, we obtain (K% = K% z $ = for I = 1, . . . , 4 and x G -D. 
This completes the proof of the claim. 

We remark that the uniqueness of the ground state of H for each hole number is not 
proved at present. We hope that this will be clarified in a future study. 

6 Estimation of a Lower Bound for fig 

In this section we estimate a lower bound for We will show later that 

— > 7a,v, (6.1) 
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with 7a,tv 




, for N > (8|D|)/9. It follows from this inequality that 



I^A,N > a lA,N 



N p + 1 Ex,* (%, bl,ab x>a %) 

2\D\ 2 ($p,$ P ) 



(6.2) 



Here, we have that 



E E ( $ P > & i,A,- $ p) = E E ($ P ,e&(A04,A,<A>)> 



(6.3) 




and it is easy to find that the right-hand-side is bounded from below by 



N/2 



sE^X^p), 



(6.4) 



i=i 



where eb(k^) < e^A^ 2 )) . . . < €b(k^ D ^) is an arrangement of eb(k) with k 6 /C in ascending 
order. Substituting this lower bound into (16.21) and taking the limit, we obtain (I4.10p . 

In what follows we prove inequality (16. ip . By the spin-rotation symmetry for $ p , it 
suffices to consider the case of a =f- By the translation symmetry, we can also assume 
x e D c without loss of generality . We first show that the left-hand-side of (16. ip with a =| 
is rewritten as 

J2scD c ;x<£S,\S\=N p W X (S) _ J2scD c ;x(S,\S\=N p W X (S) ^ 

2 T,SCD C ;\S\=N P W X {S) 2(J2scD c ;x<£S,\S\=N p W X (S) + J2scD c ;x€S,\S\=N p W X (S)) 

with the nonnegative weights 




(6.6) 



where (fi z = (b Z) \a Zj \ + a z ^b z ^). To see this, note that 



c = E Ki a ^ + E b z,i a zA = E ( & *,i a «,t + °*a 6 *,T) = E 



(6.7) 




Then, since (4> z ) 2 ^/q = (which follows from a z ^a z ^o = 0), <3> p is expanded as 




ScD e 
\S\=N P 




(6.8) 
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Since a x ^ (iLes^z) — for x G S (which again follows from a z ^a z ^Q = 0), and 

II 4>A a]..- <>, .■!>,. .!>,.. ( II $*) *o\ = for S> ^ S, (6.9) 



we have that 

* P ,aWO.,i*p) = w) 2 e (f n *o, < T ^,T^i^i (n <m *o) . 

scd c \ \zes J Vzes / / 

x<£S,\S\=N p 

= 1(N P \) 2 E W X (S). (6.10) 

Z ScD c 

x<£S,\S\=N p 

To get the second line we used 

(Vi*o, < T a»,tlti*o) = 1 = (V4*o, V>I*o) /2- (6.11) 

Likewise, we have that 

($ p ,6l a 6, a $ p ) = (iV p !) 2 2 (6.12) 

ScD c ;\S\=N p 

which together with 1 16. lOft leads to 1 16.51) . 

Before proceeding, we need to introduce some notation. For each z G D e , define D o z = 
{y | \y — z\ — 1, y G -D Q }, which is the collection of the nearest neighbour sites of z. We say 
that z and z' in D e are connected if D o z PI A, 2 / 7^ 0. For S C D c which does not contain x, 
we call z an isolated point in S if 2 is not connected any other sites in S U {x}, and write 
D X (S) for the collection of these isolated points in S. It is noted that, if y G D X (S' U {y}), 
the weight W^S" U {y}) is reduced as 

^(S'UW)^^'), (6.13) 

since a y , a with \y' — y\ < 1 commutes with fe^fex,; ILes' 0^ an d thus 

(tf 0, $A*o) = 7 E E (*o, aJ i(T 4',- CT S',--S,^o) = ~ (*o, *o) • (6.14) 
4 y'eD <t=t,J. z 

Is/'— 3/1=1 

(Recall 116.111) .) We denote by V x (N p ,l) the collection of subsets 5 of _D e such that x S, 
\S\ = N P and \D X {S)\ =1. 
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Since the value of \D X (S)\ is determined for each S C D c , we have 

y: w x (s)=y: e ^)>e e w x (s). ^ 

ScD c 1=0 SeV x (N p ,l) 1=1 SeT> x (N p ,l) 

x<£S,\S\=N p 

Now fix I > 1. Noting that there are / isolated sites in S G V X (N P , I), we find 

E W X (S) = j E E W x (S) x [y e D X (S)] 

S<=V x (N p ,l) SeT> x (N p ,l)y€D e 

= y E E E w.(5)xb/G^.(5)M5' = 5\{y}] 

SeT> x (N p ,l) yeD c S'£V x (Np-l,l-l) 

= ^ E ^(5") E E xbG J D,(5)]x[5 = 5'U{ 2/ }] 

Z< S'e2? !C (JVp-l ) J-l) SeV x (N p ,l)y£D e 



p V / S'ex»o:(iVp-i,z-i) 



To get the second line, note that removing an isolated point in S G T> x (N p ,l) yields an 
element in T> x (N p — 1, / — 1). The third line follows from (I6.13p . The last inequality is obtained 
as follows. Each site z in D e has 8 connected sites. Therefore, for every 5" G T> x (N p — 1, / — 1), 
there exist at least |D e | — 9N P sites, y, such that y is an isolated point in S' U {y}, and 
5" U {y} becomes an element in V X (N P , I). Note that \D e \ — 9N P is a positive number by the 
assumption. Then, by using I < N p , we get the last inequality. 
From (16.151) and (I6.16P we get 

E W X (S) > - 9N P ) jr J2 W X (S). (6.17) 

ScD c ZiV P \ Z / 1=0 S£V x (N p -l,l) 

x<£S,\S\=N p 

Here, for x ^ S, we have 



= ^ ( (n *o, (i - wix^ (n *°) 

< (6.18) 
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It follows from this inequality and (16.171) that 



ScD e 
x<£S,\S\=N p 



E w a (S) > 




9) E w x (s). 



- 9Ar P E E w x (su{x}) 



/ ScD c 



/ 1=0 SeT> x (N p -l,l) 



(6.19) 



xeS,\S\=N p 



From (16. 5p and (16.191) we obtain the desired inequality (16.11) . 

7 Summary and Remarks 

In this paper, for the even numbers iVh of holes in \D\ < < 2\D\, we have constructed 
a pairing state $ p with ci-wave symmetry which is expanded in terms of the Zhang-Rice 
singlet states. We have calculated upper and lower bounds of the ODLRO parameter for $ p 
as a function of the hole concentration. We have also presented the concrete Hamiltonian 
H = H + H 1 (15.61) on the Cu0 2 plain which has $ p as its ground state. We have proved that 
the lowest energy states of H (15 .7p are the Zhang-Rice singlet states and then have shown 
that, by using the positive-semidefiniteness of H 1 (15. lip , the pairing state <3> p consisting of 
the Zhang-Rice singlet states attains the ground state energy of the whole Hamiltonian H. 
The uniqueness of the ground state is not proved at present, and we leave this as a problem 
in a future study. 

It is noted that Ho with Jo = becomes the Hamiltonian of the d-p (or 3-band) model 
in the atomic limit [31 HI E], and Hq with Jo 7^ is essentially the same as the effective 
Hamiltonian derived by taking into account the hopping terms between Cu- and O-sites as 
a perturbation in the limit The idea of the Zhang- Rice singlet is based on this effective 
Hamiltonian, and the t-J model is obtained by furthermore considering the motion of the 
Zhang-Rice singlets perturbatively with the inclusion of the antiferromagnetic interactions 
between Cu-holes 



which is the effective interaction due to the hopping process between neighbouring Cu- 
sites [S]. 

In the |-D|-hole case, the present Hamiltonian has degenerate paramagnetic ground states 
with one hole per Cu-site and does not exhibit antiferromagnetism which is essential to high- 
T c cuprates. This will be improved if we consider the modified Hamiltonian Hq + Hi + H2. 
This Hamiltonian or more generally the d-p Hamiltonian with Hi may be able to reproduce 



H 2 — J 2 E S x ■ Sy, 



(7.1) 



x,ydD\\x— y\=\ 
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the essential features of high-T c cuprates, such as antiferromagnetism at low doping concen- 
trations and charge density order (or a stripe structure) between the antiferromagnetic and 
the superconducting states. We believe that further investigations about modified models 
based on our Hamiltonian which is now shown to exhibit ODRLO with <i-wave symmetry 
will contribute the understanding of high-T c cuprate superconductivity. 
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A Appendix 

In this appendix we shall show that the pairing state $ p is non-vanishing when the number 
of holes, N h = \D\ + N, satisfies N = \D\ - 2l 2 L x with some integer < l 2 < (L 2 - 2)/2. A 
similar argument will show that $ p is non- vanishing for 2L\ < N < \D\. 

It is easy to see that the collection of the states in the right-hand-side of ( 16. 8 j) is orthog- 
onal. So $ p is non-vanishing if one of those terms is non-vanishing. We shall show that this 
is the case. Let 

A\ = {x = (x\, x 2 ) | 1 < xi < Li, 1 < x 2 < 2l 2 , x 2 is odd} (A.l) 

and 

A 2 = {x = (xi,x 2 ) | 1 < x\ < Li, 1 < x 2 < 2l 2 , x 2 is even} . (A. 2) 

Now we pick up the state in (I6.8P corresponding to the subset So = (Ai U A 2 ) fl D c . Substi- 
tuting (f) z = a z ^b z ^ + &z,|0«,f into this state, we obtain 



n ( a *A,t + Ki a z,d^o = 2 ( n ^a.t n 

zeS TcSo \zeT J \z£S \ 



b zd a zA ) tf - (A.3) 

The collection of the states in the right-hand-side of the above expression is again orthogonal. 
Let Si = Ai fl D c . Then it is easy to see that 

n a z,i) n q *,t ) *o, n ^aa n ^aka ^o) (a.4) 

K zeA! J \z£A 2 ) \zeSi ) \zeSo\Si J I 

is non-zero. This implies that the term in ( 1A.3I) with T = S\ (and thus the term with S = Sq 
in ( 1A.3I) ) is non- vanishing, which concludes that $ p is non- vanishing. 
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